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Abstract. The shell model Monte Carlo (SMMC) approach allows for the microscopic 
calculation of statistical and collective properties of heavy nuclei using the framework of the 
configuration-interaction shell model in very large model spaces. We present recent applications 
of the SMMC method to the calculation of state densities and their collective enhancement 
factors in rare-earth nuclei. 


1. Introduction 

Collective states constitute a significant part of the spectra of heavy nuclei up to moderate 
excitation energies and they are often described by bands that are built on top of intrinsic states. 
However, the contribution of these collective states to the nuclear state density is difficult to 
calculate microscopically, and is often described by phenomenological collective enhancement 
factors pp. The challenge of computing the collective enhancement factors microscopically can 
be overcome using the shell model Monte Carlo (SMMC) method [21 [31 SI E]- This method is 
particularly suitable for the calculation of thermal and statistical properties of nuclei. 

Here we discuss recent applications of the SMMC method to the isotopic families of even- 
even samarium and neodymium nuclei m E]: which exhibit a crossover from vibrational to 
rotational collectivity. In particular, we present results for the state densities and their collective 
enhancement factors. 

2. The shell model Monte Carlo (SMMC) approach 

The SMMC method is based on the Hubbard-Stratonovich (HS) transformation |8] to express 
the Gibbs operator of a nucleus (described by a Hamiltonian H at inverse temperature 

j3 = 1/T) as a superposition of one-body propagators of non-interacting nucleons moving in 
external auxiliary fields <T(r) that depend on imaginary time r (0 < r < /3) 

= I D[a]G^U„ . ( 1 ) 

Here is a Gaussian factor and Ua- describes a one-body propagator associated with a given 
set of auxiliary fields a. Subsequently, the thermal expectation value of an observable O at 


inverse temperature (5 can be written in the HS representation as 


^ Tr(Oe-/^^) ^ jD[a]WMOU 

\ / Tr{e-f^H) jD[a]W„^a ’ ^ ’ 

where (O)^ = Th: {OUcj)/ ^ is the thermal expectation value of the observable in a given 
configuration of the auxiliary fields a. Since the numbers of neutrons and protons are fixed for a 
given nucleus, all traces in Eq. ([2]) are evaluated in the canonical ensemble. Dehning a positive- 
definite function VFo- = Go-|Tr Ua\ and the associated Monte Carlo sign <ho- = Tr Ucr/\Ti C/o-|, 
auxiliary-held conhgurations ak are sampled according to Wa-, and the expectation value in ([ 2 ]) 
is then estimated from (O) « Y.k{0)at,^ak/Y.k 

3. Collective enhancement in the state densities of heavy nnclei 

Nuclear collectivity is a ubiquitous phenomenon in heavy nuclei. Various types of collectivity, 
such as vibrational collectivity and rotational collectivity, are observed and are well described 
by phenomenological models. However, a microscopic description of nuclear collectivity using 
the conhguration-interaction (Cl) shell model is a challenging task. In particular. Cl shell model 
calculations of heavy nuclei require model spaces that are many orders of magnitude larger than 
model spaces that can be treated by conventional diagonalization methods. The SMMC was 
successfully applied to describe the rotational character of the rare-earth nucleus ^®^Dy [9] in 
the framework of a truncated spherical shell model space. Recently, using the same approach 
and model space, the crossover from vibrational to rotational collectivity in families of even- 
even samarium and neodymium isotopes was also reproduced Here we discuss the state 

densities of these nuclei and the associated collective enhancement factors. 

The single-particle model space consists of the 50 — 82 major shell orbitals plus the I/ 7/2 
orbital for protons, and the 80 — 126 major shell orbitals plus the Ohii /2 and lgg /2 orbitals 
for neutrons. The bare single-particle energies of the orbitals were chosen so they reproduce 
the Woods-Saxon energies in the spherical Hartree-Fock approximation. The effective two- 
body interaction consists of monopole pairing interaction terms for protons and neutrons, 
and multipole-multipole interaction with quadrupole, octupole and hexadecupole terms. The 
calculations were performed using an SMMC code in the proton-neutron formalism [9]. 

3.1. State Densities 

The SMMC method has proven to be a powerful technique for the calculation of state 
densities mm- In the SMMC approach, the thermal energy is calculated as the expectation 
value of the many-body Hamiltonian at an inverse temperature (3. The canonical partition 
function is computed by integrating the thermal energy with respect to /3, and the state density 
is obtained by inverting the Laplace transform of the canonical partition function in the saddle- 
point approximation. 

In Fig. [1] we show the SMMC state densities calculated for the even-even and 

144-152]^^^ isotopes. We compare the SMMC results (circles) with the experimental state 
densities at low energies using level counting data (histograms), and at the neutron resonance 
energy (when available) assuming the spin cutoff model with rigid-body moment of inertia 
(triangles). We also show the BBF state densities (solid lines) whose parameters (i.e., the single¬ 
particle level density parameter and the backshift parameter) are determined from level counting 
and neutron resonance data. The SMMC results are also compared with the densities calculated 
by the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation (dashed lines). The 
observed enhancement of the SMMC densities over the HFB densities can be attributed to the 
presence of collectivity in these nuclei, since collective excitations are included in the SMMC 
approach but are absent in the HFB densities. The HFB densities also display “kinks” at 
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Figure 1. State densities in the even-even and isotopes. The SMMC 

densities (circles) are compared with level counting data (histograms) and neutron resonance 
data (triangles) when available. Also shown are the BBF state densities (solid lines), which 
are parametrized using the experimental data, and the HFB densities (dashed lines). The thin 
arrows indicate the neutron and proton pairing phase transitions, and the thick arrows indicate 
the shape phase transitions. Adapted from Refs. mm- 


energies that correspond to the phase transitions associated with the breaking of the proton and 
neutron pairs (thin arrows) and with the transitions from deformed to spherical shapes (thick 
arrows). These shape transitions occur in all nuclei except ^^®Sm, ^^“^Nd and ^^®Nd, which are 
spherical in their ground states. The shape transitions in ^®^Sm, ^®^Sm, ^®^Nd and ^^^Nd occur 
at excitation energies that are higher than those displayed in Fig. [T] They can be seen, however, 
in Fig[2l 

3.2. Collective Enhancement 

The microscopic calculation of collective effects in the state density of heavy nuclei has been 
a major challenge. Such effects are usually described in terms of phenomenological vibrational 
and rotational collective enhancement factors [T]. 

Recently we introduced PE! the ratio of the SMMC and the HFB densities, i.e., K = 
/Osmmc/phfb, as a microscopic measure of the collective enhancement factor. In Fig. O 
we show R as a function of the excitation energy for the above families of samarium 
and neodymium isotopes. In the spherical nuclei ^^®Sm, ^^^Nd, and ^^®Nd, collectivity is 
predominantly vibrational and is observed to be lost completely (i.e., AT ~ 1) above the pairing 
transition energies. The other samarium and neodymium isotopes are deformed, and collectivity 
persists at higher excitation energies than the pairing transition energies where the enhancement 
is due to rotational collectivity. Above the shape transition energy, the collective enhancement 
factor decays to AT ~ 1 since spherical nuclei can no longer support rotational bands. 

4. Conclusion 

We discussed recent SMMC applications to the even-even and isotopes. 

The total SMMC state densities were shown to be in very good agreement with experimental 
data. We also extracted a microscopic measure of the collective enhancement factor dehned 
by the ratio of the SMMC and HFB densities. We observed that the decay of vibrational 


















100 


\ ■' 

■ 

'^®0 0 o n 

h 1 

‘p ° ° ° < 

j^®o®ooj . 

|n 0 


148 

Sm 

150 

Sm 

152 

Sm 

154 

Sm 







ir 



O 

© 

o 

o 

0 

0 

0 

® 0 

144 

Nd 

<*hH¥a«^<3BOQ0g O O 

146... 

Nd 

148.TJ 

Nd 

150.. . 
Nd 

152.. . 
Nd 


0 10 20 0 10 20 0 10 20 0 10 20 0 10 20 30 

(MeV) 


Figure 2. Collective enhancement factor K versus excitation energy Ex in the even-even 
148-I54gj^ and isotopes. The pairing and shape transition energies are shown by 

the thin and thick arrows, respectively. Adapted from Refs. [61E!. 


and rotational collectivity with excitation energy correlates with the pairing and shape phase 
transitions, respectively. 
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